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INTRODUCTION 

This is the first quarterly report on a one year program to 
provide computational methods for prediction of solar cell perfor- 
mance in a natural radiation environment. It covers work performed 
during the period 1 June 1969, through 31 August 1969. A model, 
permitting theoretical calculations for solar cell performance 
under such a mixed-radiation environment, was previously prepared 
under contract 952246 with the Jet Propulsion Laboratory, Pasadena, 
California. That model, reported by Exotech on February 28, 1969, 
»is the basis for the effort reported herein. The successful 
adaptation of such a model in a computerized version sensitive to 
environmental parameters is clearly an invaluable aid in design, 
selection, and performance predictions for individual spacecraft 
power systems. 

For readers unfamiliar with the previous work. Section I 
provides a brief synopsis and review of the mathematical terminology 
introduced. Section II and III are really the substance of this 
report;- they contain the progress achieved in the past three months. 
It will be obvious, on reading these sections, that portions of the 
model have been adapted to a computer to permit the investigations 
reported. To Indicate the extent of this adaptation, Section IV 
presents the current status of the program. Finally, conclusions 
regarding this work and observations of potentially fruitful areas 
for photovoltaic research are presented as Section V. 
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I. SUMMARY OF SOLAR CELL PERFORMANCE 

The photovoltaic effect at a silicon diode junction was dis- 
covered as early as 1941, ^ but practical- silicon solar cells were 
not available until 1954. These solar cells converted sunlight with 
an efficiency around 67„ and opened up the possibility of rugged, 
long-lived, light power sources for special applications. An interest- 
ing review of early developments of the solar cell is given by 
Crossley, Noel, and Wolf . ^ 

Modern solar cells for spacecraft applications are typified by 
Figure 1. The cell area may be from 1 x 2 cm to 3 x 3 cm, with a cell 
thickness of 8-12 mils (dimensions are conventionally given in these 
heterogenous units). A junction is formed about 0.5 micron (5 x 10~ 5 cm) 
below the sunward surface; this junction separates the thin "surface" 
region of n-type silicon from the thicker "base" region of p-type 
silicon. For the sunward side, the electrical contact is in the form 
of a bar along one edge and a number of thin grid lines extending from 
it across the surface. The dark side contact generally completely 
masks the back surface of the cell. 
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■ Surface region (about 0.3 p,) 
1* — Base region (about .03 cm) 


-> x 


Si 


Figure 1. Exploded view of typical solar cell, with notation 
as used in text. 
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Illumination of this device results in the sunward surface 
being a fraction of a volt positive with respect to the back surface. 
A current is thereby induced in an electrical connection between 
the front bar and the back contact. This is not equal, however, to 
the current generated by the sun (photovoltaic current) , for some of 
the photovoltaic current is returned through the solar cell itself, 
in accordance with its diode property. This leads to the solar cell 
equation, which principally states that the electrical current flow- 
ing from a solar cell equals the difference between the photovoltaic 
.current I produced in it, and' the diode current T lost in it. 

i " D 

The output of the solar cell depends on a large number of in- 
dependent factors. First, the intensity of sunlight and the angle of 
incidence are important. A cell directly facing the sun receives the 
maximum possible sunlight. If it deviates from this direction by an 
angle 0, as shown in Figure 1, the radiant energy striking the cell is 
.reduced, by the cosine of 0. As the angle increases, edge effects, 
especially when a coverslide is used, make the reduction deviate 
slightly from this law, and measurement must be relied on for each 
specific geometry when extreme accuracy is required for solar cells 
illuminated at large angles from the perpendicular. 

The output of solar cells has frequently been characterized by 

the short-circuit current I , the open circuit voltage V and the 

sc oc 

power and voltage of the cell near its maximum power point. While 
these parameters do characterize the electrical output, in the pre- 
sent work more basic terms are used. These are "the photovoltaic 
current I , the diode parameters I and V , and the internal series 
resistance R. These four parameters are chosen since it is a more 
straight -forward exercise to determine from semiconductor theory how 
they are affected by the environment and history of a solar cell. 
Neverless, both set of parameters are related, and knowing one set 
permits a calculation of the other. 

In this connection, it is helpful to note that I nearly 

sc 

equals I , that R is generally less than 0.5 ohms, and that V is on 
L o 


the order of 40 millivolts. An example of numerical values in 
the solar cell equation is presented in section II. 

Aside from the obvious factor of illumination, the other 
factors affecting solar cell performance include its temperature T 
and the fluence of atomic particles $ to which it has been exposed. 

The fluence of natural radiation in the space environment is 
composed of many species of atomic particles, but only protons and 
electrons ' need to be considered for their damaging effects. To 
determine the effects of this radiation on solar cells, it is 
’necessary to (a) evaluate the proton and electron fluences as funct- 
ions § (E) and $ (E) of energy, (b) evaluate and sum the damage 

P £ 

due to these fluences that occurs in the solar cell material, (c) 
evaluate the changes in the solar cell electrical characteristics as 
functions of this damage. 

The current I at a voltage V across the solar cell may be 
evaluated from its radiation history, as described above, and a 
knowledge of its construction, as indicated by Figure 1, the illumi- 
nation,' and the temperature. This is the task of our mathematical 
model. 
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. II. ANALYSIS OF MATHEMATICAL MODEL 

A. Comparison with Recent Literature 

A continuing search of the current literature is being 
undertaken in order to determine the existence of possible variations 
to the mathematical model we have adopted ^ for solar cell perfor- 
mance. This search indicates that validity of the diode equation for 
silicon solar cell applications is being tested by several investigators. 

The current-voltage characteristic of the solar cell is assumed 
„to be equivalent to that of a current source, shunted by a diode, and 
with a series resistance. The equation resulting from this assumpt- • 
fon is the conventional solar cell equation. 1 * We have developed 
curve-fitting techniques that so far have shown this equation to fit 
experimental curves with a high accuracy under all conditions of 
illumination and radiation degradation. Similar sucess has been 
reported by Brown. ^ ^ To demonstrate the accuracy possible with the 
simple solar cell equation even for unconventional cell types, we 
present two curve fits in Figure 2 to lithium-doped solar cells. 
Generally, our curve fits have been almost, but not quite, this 
good. 


Theoretically, the solar cell equation appears to be a good 
first approximation. This is borne out by its accuracy in fitting 
current-voltage *(1V) measurements. However, deviations of measure- 
ments from the equation have been reported. Wolf and Rauschenbach^ 
have recommended that the current source be assumed to be shunted-by 
two diodes with different characteristic voltages. The revised 
solar cell equation, in our notation, would then take the form 



(V+IR sl> /V ol 


ol 


- 1 i 


(V+IR s2>/ V o2 


o2 



( 1 ) 


This has the immediate advantage of allowing two extra 
parameters for curve fitting, and the disadvantage of relating these 
parameters to the individual solar cell and its environment. The 
rationale for the extra diode term must be developed if it is to be 
adopted in a general mathematical model. 
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V (millivolts) 


Figure 2. 

Comparison of measured IV curves for 
with points computed with the solar 

lithium cells 
cell equation: 

I = 60.2 - 

1.675 x 10 -V''-K791I)/45.34. 1) 

top curve 

I - 46.3 - 

1.625 x 10- -* (e <V+.322DM6,»8 

bottom curve 






Schoffer and Beckman studying the response of solar cells 

to intense illumination, derived this model from theory. Much of the 

series resistance R in the cell is due to the sheet resistance of the 

s 

surface layer, through which the current must flow from all points of 
the junction area to reach to grid lines. With large currents, this 
results in a appreciable voltage gradient on the junction area, with 
IR drops in voltage between points at the grid lines and points far 
from them. Thus, a distributed diode potential exists, and the two- 
diode model of Wolf and Rauschenbach can be considered as an approxi- 
mation to the distribution. Further, measurements ol open circuit 
voltage as a function of illumination plot as a slightly curved line 
in a semilog plot; instead of a straight line as predicted by the 
simple solar cell equation. This is interpreted as indication that, 
with a change in the IR drop along the surface, the voltage across 
the diode changes, and the relative magnitudes of the two diode cur- 
rents shift. 


This last point is worthy of a more detailed discussion. 

Consider a solar cell with negligible resistance R , which obeys 

* ■ s 

the simple solar cell equation. Then 1^ equals l gc > and 


I - X -I ( e V/V °-X ) 
sc o v 7 


( 2 ) 


This rearranges to 


I - I » 
sc 


v 




(3) 


Since V q is typically about 50 millivolts, when the voltage 
V is greater than 350 millivolts, e V ^ /o is greater than 1096 and, to 
better than 99 . 97 % accuracy. 



I e 
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V/Vo 
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. For such a solar cell. This means that a plot of the logarithm 
of ^sc" 1 ^ versus the voltage V, beyond the knee of the I-V curve 
which generally occurs around 350 mv, should be a straight line 




Figure 3, A plot of the diode current (I - I) from data taken 

/ON SC 

with different light intensities. ’ The plots are straight 


where the data obey a simple, resistance-less solar cell equation. 
Resistance would make the lines concave upward. The observed 
downward concavity suggest the existence of additional diode effects 




with slope determined by V q and extrapolated intercept determined 

by I . 

° 

An attempt at such a result is shown in Figure 3. The lines, 

drawn for different illumination intensities, are almost straight, 

but not quite. It appears that when the diode current (I ~ I) is 

s c 

less than about 40 rav, the lines are parallel and straight* This 
has the effect of making some solar cells deviate from simple diode 
theory when the illumination is great, as was discussed. 

Whatever the reason, this indicates that the simple solar 
« cell equation is only an approximation. The model is inaccurate 
when the diode current is large; 40 ma in a 2 cm cell appears to be 
the borderline. 

It is not necessary to invoke distributed diodes to derive 
a curvature such as exhibited in Figure 3. Ladany^ 9) has presented 
an analysis of a one-dimensional diode, such as the solar cell, 
which results in a forward characteristic reminiscent of that shown 
3. His analysis leads to an expression, however, that 
does not require additional parameters to duplicate the characteristic. 

B. Computation Techniques 

The heart of the Kxotech mathematical model is a difference 
solution for the continuity equation. This allows a calculation of 
the minority carrier profile across the cell for arbitrary light 
spectrum and for nonuniform damage by radiation. 

The difference technique employs a mesh interval h whose width 
determines the convergence of the solution. It appears reasonable 
to assume that the smaller h is near the junction, the better-defined 
is the profile there. Since the minority carrier distribution near 
the junction determines the photovoltaic current, it is desireable " 
to calculate it as precisely as possible. To demonstrate this 
convergence, we present Figure 4 from earlier work. The figure 
suggests that decreasing the mesh interval improves accuracy, to a 
point where h is about a quarter of a micron. Whether there is an 
optimum mesh interval, as suggested by the figure, must be .determined 
by more detailed comparisons. 
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h (microns) 


Figure 4. Calculated values of short-circuit current density versus 

mesh interval h, for three cells with different base minority 
■ ’ carrier diffusion lengths L. 

Small mesh intervals such as these are undesirable in terms of 
computer time to solve the continuity equation./ If h is uniformly 
0.25 microns, crossing the base region of a typical 10 rail cell would 
require about 1000 intervals. Each iteration of the calculation would 
have that many steps. . 

However, such precision is not necessary near the back of the 

solar cell. The ideal scheme for small mesh size near the junction, 

yet reasonably rapid computer time, would therefore be use of a 

formula by which a mesh interval h increases with distance x from-" 

K k 

the junction. We selected the formula 

h k = kA (5) 

If, then, the distance b is to be divided into n Intervals, the first 
interval will have a width 


10 , 




A = 2b/n (n + 1) 


( 6 ) 


and the last or nth interval will have width nA. As an example, with 
only 50 intervals in a 10 mil distance, the first has a width 0.2 
microns. 

To utilize this scheme, it is necessary to rederive the 
difference equation for a mesh of variable width h^.. The first and 
second derivatives of the minority carrier concentration, designated '<■ 
n^ at the beginning of the kth interval, must be evaluated. The 
approximate derivation we used is based on Figure 5. 



Figure 5, Construction of variable mesh hu to approximate 
the curve n(x) at points n^ 

Figure 5. Construction of variable mesh h to approximate the curve 
, K 

n(x) at points n^. 

Approximating the curve n(x) by straight line segments joining 
points n^. and using the definitions of the derivative in calculus, 
we have 

( f! ) ‘ [ Vl + ( ) ( "k+2 - ”W - "k ] n \ (7) 

b h k+l ' 
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where b denotes where the first derivative is computed: this is at 
m k+l . 


Precise calculation of the second derivative at n^_^^ requires 
more points than the three in the figure. This leads to difficulties 
in the continuity equation. However, it can recalled that h^ is 
generally small, so that the first derivative changes little within 
each interval. From to in the figure, the first derivative 

is approximately 

(£) - 3 .<« 

and from n^ +1 to n^ the first derivative is approximately 



= n k+2 " n k+l 

Vi 


m 


So that the second derivative at 


it j_, is approximately 



bfc "k+2 Wl* Vl + Vl ”k (10) 
2h ^ Vl 


These expressions for the first end second derivet ives trans- 
form the continuity equetion into e differense equation, with vari- 
able mesh size. After some routine shuffling, the terms can be 
organized to yield 


V2 


\ + h k j “ Vu [ \ ( 1 +h k i E /«) + Vi d-tyjE/kT) 


[ Vl O-VE/kT -2h= /l£ ) 
" L 2h k Vl G k-H /D ] 


( 11 ) 



This replaces the earlier formulation of the difference equation 
to permit faster, more accurate calculation. 

C. Electron Shielding Calculation 

An elaborate curve fit was presented in the mathematical 
model to relate the spectrum of electrons in the natural environment 
with the spectrum striking a solar cell that is covered by a fused 
silica coverglass. This fit has been programmed in the last quarter 
and accurate calculation of the transmitted spectrum was found to 
require an undesirably long and detailed calculation. 

The difficulty is in the shape of the spectrum formed when 
monoenergetic electrons are slowed down in a coverglass. A typical 
spectrum, as seen in Figure 6, includes most of the transmitted 
electrons in an almost monoenergetic peak. This peak decreases, 
and the transmitted spectrum becomes flatter as the coverglass thick- 
ness approaches the range of the incident electrons. 

We found that the curve fit to the spectrum can not be integrated 
by ordinary analytic techniques. This forced us to program the 
computer to calculate points on the curve, and generate damage integrals 
by a piecewise approximation 

<K * > = X W < E 1 ’ V < 12 > 

k • \ 

where (E ) is the electron damage coefficient for transmitted 

energy and (E , x fc ) is the flux, evaluated from the curve 

fitted spectrum, of energy due to monoenergetic electrons of 

energy E slowing down through a shield of thickness x, . (The 
o k 

notation here is unchanged from the previous report.) 

Strictly speaking, $ e (E^ , x k ) should be a piecewise integration, 
covering all of the spectrum in an energy interval about E.. 
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If the spectrum were reasonably flat, this could be approximated by 
calculating the energy-differential flux at energy and multiplying 
it by the width of the interval. The peak in the spectrum makes this 
an impractical approach. . , 

Our most recent programming has followed an idea previously 
presented. The redefined formula is 

r ' ‘ 

< K:{> k =£. K e (E t . x K ) i e <E t ) ( 13 ) 

where is now an energy of the incident, or natural environment, 

electron spectrum and the electron damage coefficients are modified by 

the shield thickness x, . 

k 

2 

When the solar cell is covered with pd grams/cm of coverslide, 
the electron flux due to incident monoenergetic electrons of energy E is 
given by a spectrum, of energies less than E, dependent on the thickness 
pd of shielding. This spectrum can then be weighted by the damage 
coefficient as a function of energy to obtain an effective damage 
coefficient K e (E , pd) for the shielded solar cell with electrons of 
energy E incident on the shield. 


Several results of this calculation are shown in Figure 7- where 


Figure 7. Effective damage 
coefficient of electrons 
as a function of coverslide 
thickness (gms/cm 2 ) and 
incident energy 
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the damage coefficient, per unit incident particle, appears to decrease 
exponentially (e~ lipd ) . When pd approaches the electron range, of course, 
the proportion rapidly drops to zero. Further, the fitting parameter p> 
appears to have a simple dependence on energy. The slopes of the 
curves are plotted, in Figure 8. versus incident particle energy. 



Figure 8 . Plot of negative slopes of straight lines in Figure 7 versus 

electron energy, compared with negative slopes that would yield data of 

2 

Ref. 12 for solar cells shielded with a 0.3 gm/cm coverslide (See Eq. A) 



The results fit well to an analytical expression: 

K e (E,pd) = K (E ,0) exp (-lOpd/E 1 / 5 ) (14) 

' 2 
where pd is the coverslide areal density, in gm/cm . The equation 

14 is obviously adaptable to the problem of determining the necessary 

coverslide thickness for solar cells on a given mission. K is now 

the effective damage coefficient, having been reduced in magnitude 

by the presence of the coverslide. 

We have compared our results with the early experimental work 

( 12 ) • 

e of Brown et al Using a beam at different angles on bare and 

variously-covered cells, they constructed the damage curves for n/p 

solar cells for electron energies up to 3 MeV, and extrapolated to 

7 MeV. Assuming that the damage is given by the data for bare cells 

2 

and the data for cells covered with 0.3 gm/cm coverslides, and by 
a simple exponential exp (-ppd) at other thicknesses, we have computed 
Ms a function of energy. The energy dependence of the function p. 
is also shown in Figure 8. The dependence on E again indicates a 
minus 1.5 power law; the magnitude of p- is however, about 30% greater. 
Part of this discrepancy may be due to the decrease in ionization 
loss by electrons with higher atomic weight of the shield. 



III. PROTON DAMAGE STUDIES 

A. Depth Dependence of K 

Because' a proton slows down as it travels to a stop in a 
medium, its damage effectiveness increases. The result is that K may be 
considered as a function of depth x into a solar cell such that the 
change in diffusion length L varies with depth. The angular distribut- 
ion about the normal to the surface describing the incident proton flux 
affects this distribution. As an example of the magnitude of the 
variation with depth, consider the case where the incident flux is 

isotropic and the. variation of K with energy E is simple K /E where K 

o o 

is the damage coefficient for 1 MeV protons. 

Let equal the cosine of angle ® of incidence, as shown in 
Figure 9. The energy of a proton in its travel can be determined from 
the range R^E^ of the proton. By integration over all possible 

angles of incidence, the express- 
ion for K is obtained. 



K(x) = 2u 


K 


Mu 


4rr 


E n - 
. o 


X 

R p* 

o 


n 


1 dp, (15) 


Figure 9. Proton stopping in 
solar cell 


where is the cosine of the maximum 
angle of incidence through which a 
proton of energy can penetrate to 


the depth x. 

The expression 2 tt represents integration over the azimuthal angle, and 
the incident proton flux is normalized to (1/4tt) protons per steradian. 
With the substitution 


E 


n 


x 

R p 
o 


(16) 


the integral is transformed to 

K (x) 


(E n - x/R ) 

p o o 


K(x) = 


2R 


J 


dz 


= 1/n [E n - z] 2 


(17) 
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Even with the simplifications presented so far, this integral is 
formidable. With n equal to about 1.75, 1/n equals 0.57. As an 
approximation, let 1/n equal 0.5. Then a closed form solution is 
obtained: 


K 

K(x) - 

o 


(1 - 2 ) % 
R 



( 18 ) 


where R is the initial range of the protons. 

«, The various parts of the expression may be interpreted. The 

expression K q /E o represents the damage coefficient of the incident pro- 
tons. The factor 1/2 represents the shielding of the back hemisphere 
by the cell. The expression in brackets represents both the reduction 
in number of the protons with depth x and the increasing damage of the 
remaining protons as their energy drops. This expression is plotted in 
Ris^^® 10. It can be seen that the drop in energy at first enhances the 
damage, but the drop in number eventually takes over and reduces the 
damage until it vanishes at the point where depth x equals proton range R. 
This is quite different from the damage profile of monoenergetic protons 
entering the solar cell in a parallel beam. The damage is "smeared" 
for the isotropic incidence; it is "peaked", as suggested by the • 
familiar Bragg curve, for the parallel beam. 

B . Model Predictions 

A major outstanding problem in the prediction of radiation 
effects to solar cells is the calculation of damage from protons that 
do not completely penetrate the crystal. . The model developed under 
Contract 952246 provided an original method for calculation of the 
effects of nonuniform damage in solar cells, such as generated by low 
energy protons. Thus, these proton damage relationships could be 
evaluated for the first time in depth. Previous work has been based 
on "equivalance" or "effective" measured parameters. 
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Figurell. The damage due to unit isotropic fluence, normalized 
to the surface damage due to a unit parallel-beam fluence of protons. 
Plotted. as a function of depth of penetration into the solar cell, with 
the depth normalized to the proton range (Approximate, from Eq. 18) 
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As a first possibility, it might be possible that our value of 
the damage coefficient is too high. The measurement of proton damage 
coefficients is obviously difficult, in view of the rapid change of 
proton energy even in s small silicon crystal. Further, the coe- 
fficient will depend on crystal resistivity, type of dopant, and even 
temperature during bombardment. These dependences are poorly understood 
and no simple formulas exist to represent them. Thus, K could easily 
•have been overestimated. If,, for example, it were a factor of 3 lower 
than assumed, then the model would estimate a given decrease in I 

L 

after one- third of the proton fluence actually required. 

It is a simple matter to check such a possibility. Figure 11 
includes a dashed line for the calculation where K is reduced, for all 
proton energies, by a factor of 3. The alternate model generated in 
this way approaches the measurement, but some further scaling of K will 
show that a reasonably accurate fit is not possible from simple scaling. 

A second possibility that might be considered here is interference 
between defects ^at the end of the track of the proton. Each proton 
produces about 5 x 10 5 displacements/cm as it is slowed down from an 
energy of 10 keV to 1 keV. Conceivably, in a large fluence of protons, 

i 

many of. the collisions that occur could be with silicon atoms that 
have already been displaced. These can not be counted at additional 
defects and therefore should not add to the damage. However, under 
the assumption that all silicon atoms, whether displaced or in lattice 
positions, are equally likely to.be struck by a proton, we calculate 
that this interference should be significant when the proton fluence is 
of the order of 10^/cm^. 
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As was shown in the last section, the proton damage 
coefficient can be computed as a function of depth into the solar 
cell. The generalization of this calculation to the proton spect- 
rum in space is obvious. When a covers lide is present, the para- 
meter x in Eq. 17 must be measured from the front surface of the 
coverslide. Taking these considerations into account in the model, 
we compute the damage, and hence the values of the minority 
carrier coefficient as a function of depth in the cell. Eq. 11 
can then be solved for the minority carrier concentration, and 
'hence the current. In theory, this precise analysis is preferable 
over the practice of estimating an effective damage coefficient and 
an effective L for the cell. 

With this detailed analysis, we still find in some cases that 
theory overestimates the reduction in current. The experiment 
compared with is an extreme one: monoenergetic protons were directed 
perpendicular to a solar cell. Thus, the damage is strongly peaked 
near the end of the range of the protons, and Lj, has a strong dis- 
continuity. The depth in the cell at which this occurs in about 3.3 
microns. The results we have obtained are shown in Figure 11. 

Several possibilities exist to explain why low energy protons are not 
damaging as- predicted. 



Figure 11. Photovoltaic current versus 270 keV proton fluence. 
The dashed line is based on the assumption that the damage coefficient 
used is a factor of 3 too large. 




The deviation between theory and measurement occurs much lower, 

12 2 

around fluences of the order 10 /cm . However, the assumption made 
may be in error; the displaced silicon atoms could be more likely to 
be struck than are the lattice atoms. This would be a channeling 
effect, whereby lattice atoms, particularly those at the end of the 
proton track, shadow each other to reduce a collision probability. 

This needs a more detailed theoretical investigation. 

. Another possibility is that displacements occur and interfere 
with each other. A lattice vacancy, according to theory, diffuses 
through the lattice until it finds a stable configuration. One such 
configuration is a recombination center, which is electrically active 
in reducing the carrier density. Another configuration, of course, 
would be a vacancy filled with a silicon that had been interstitial. 
This would be electrically inactive, and be readily evaluated. in 
terms of theory. 

. '? 

A fourth possibility for investigation takes into account the 
proton, or hydrogen ion, that is left at the end of the track. The 
proton could combine either with a nearby vacancy or with a recombinat- 
ion. center to reduce its damage effectiveness. This would be local 
form of annealing, perhaps similar to that noted with lithium. 

(These atoms are similar, having the same valence, and similar size.) 

Finally, the recombination center may persist but be incorrectly 
evaluated. The current study does not include evaluation of fill 
factors, which are known to be significantly affected by , illumination 
intensity. This topic has recently been of concern and is the subject 
of current investigations elsewhere. Fill factor for proton- induced 
recombination centers is a function of light intensity, and therefore 
can vary through the cell. This is because light is strongly attenuated 
in silicon. The light intensity at the recombination center is there- 
fore a function of its distance from the surface of the cell. 

The majority of these possible explanations invoke some sort of 
interference between defects. If one of these is correct, a laboratory 
experiment with monoenergetic protons in a parallel beam should create 
less damage than the same fluence in space. As shown in the previous 
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section, the assumption of an isotropic space fluence leads to a 
damage constant K that varies smoothly through the cell, and there is 
no region of dense damage. This would decrease the amount of defect 
interference, displacement interference, or localized . annealing. As 
a result, damage would be distributed more evenly through the solar 
cell and possibly more accurately predicted by theory., 

An experiment that would shed light on the true nature of proton 

damage could be readily performed in the laboratory. This would 

consist of measurements of proton damage to solar cells that initially 

: were matched as to crystal orientation, junction depth, efficiency, 

etc. A proton energy, in the neighborhood of 1 MeV, would be selected t 

and each cell would be bombarded at a different angle with respect to 
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the beam, with a fluence on the order of 10 /cm . The current output 
of the damaged cells would then be measured under a solar simulator. 

The results would provide a test of common assumptions that the 
damage resides at the origin of the defect, that channeling is of 
no concern, etc. A careful comparison of the results with the 
model should permit a better understanding of the effects of non- 
uniform damage: and the diffusion of charges in silicon. 
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IV, STATUS OF COMPUTER PROGRAM 

Portions of the mathematical model have been programmed and 
run with test data on the IBM 360, These runs have led to the several 
modifications for the model that are discussed earlier in this report. 
Thus, the computer program at present exists as a series of routines 
subject to revision of theory and technique. Future effort is necess- 
ary to compact and join these routines. 

The central routine developed so far is the calculation of the 
photovoltaic current, via an analysis of the continuity equation. 

The variable mesh size h^ which is discussed in section IIB is now 
Incorporated in this program. The operator enters the program in the 
machine, together with data on cell thickness, junction depth, etc. 

The program can calculate minority carrier concentration on both sides 
of the junction and sum the photovoltaic currents from both sides. 

Three routines have been written to determine the damage integral 
due to trapped protons, trapped electrons, and flare protons. In the 
first tvjo of these, a fine mesh has been set for the spectrum so that 
reasonable accuracy appears possible when the spectrum is assumed 
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flat in the interval. For flare protons, a spectrum varying as E is 
assumed and the portion of the damage integral due to the high energy 
tail is performed analytically. 

Finally, a short routine has been written to compute the source 
term for minority carriers due to the absorption of light. Values 
of the light absorption coefficient of silicon are stored in this 
routine. The spectrum of space sunlight is also stored; a future 
option that would be useful would be to include other useful spectra, 
such as tungsten arc, etc. The routine presently does not include 
consideration of reflection losses. 


V. CONCLUSIONS 

A. General Summary 

Results of the first three months of an effort to prepare a 
computer program for solar cell performance in space are reported. 
Improvements of a mathematical nature have been noted: the application 
of a variable mesh interval in the difference equation, and the re- 
grouping of factors in the electron shielding calculation to avoid a 
rap idly -varying function. 

=> , Analytically, an investigation into the diode property of the 
solar cell has been conducted. This has lead to an supper bound in 
diode current, beyond which a more detailed model is necessary. 

Problems encountered are discussed. These are principally with 
regard to the monuniform damage by protons. An approximate analysis of 
the damage due to monoenergetic protons, incident isotropically, is 
presented. Arguments are developed from this analysis to indicate that 
the discrepancy between model results and laboratory experiments 
should be greater than between model results and performance in space. 

The computer program, existing as a series of independent 
routines, is in the process of formulation into a single package. 

B. Future Work 


Future work will include further analysis of solar cell 
performance in light of current theory and observation, preparation 
of the completed computer program, debugging and checkout, and 
evaluation of the program. 

The analysis of solar cell performance at present is concentrated 

on the problem- of nonuniform damage. It should be noted that only 
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when the fluence reaches the rather high value of 10 does the model 
deviate appreciably from experimental data and, furthermore, reasons 
exist to indicate this deviation may be smaller for space fluences 
than for laboratory fluences. Further effort on this problem should 
permit calculation of damage with higher values of proton fluence. 
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The completed computer program shall be capable of predicting 
the power degradation of solar cells in a combined environment of low 
and high energy protons and electrons and provide the necessary frame- 
work to optimize the solar cell coverglass combination for a parti- 
cular radiation environment with respect to cell base resistivity, cell 
thickness, coverslide thickness, and cell substrate type and thickness. 
The program is to be presented in operable form, listed in Fortran IV, 
suitable for input to an IBM 360-50. 

Standard modular debugging will be employed for program checkout. 
It will be tested and evaluated against published experimental results 
and continually updated throughout the course of this contract as 
additional experimental data is analyzed. 

C. New Technology 

After a diligent review of the work performed under this 
contract, it was determined that no new innovation, discovery, 
improvement or invention was developed. 
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